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The non-perturbative gluon and ghost propagators in Landau gauge QCD are obtained using the 
Schwinger-Dyson equation approach. The propagator equations are solved in Euclidean space 
using Landau gauge with a range of vertex inputs. Initially we solve for the ghost alone, using a 
model gluon input, which leads us to favour a finite ghost dressing in the non-perturbative region. 
In order to then solve the gluon and ghost equations simultaneously, we find that non-trivial 
vertices are required, particularly for the gluon propagator in the small momentum limit. We focus 
on the properties of a number vertices and how these differences influence the final solutions. The 
self-consistent solutions we obtain are all qualitatively similar and contain a masslike term in the 
gluon propagator dressing in agreement with related studies, supporting the long-held proposal of 
Cornwall. 
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1. Introduction 

The non-perturbative gluon and ghost propagators may be obtained from their Schwinger- 
Dyson Equations (SDEs). They are interesting quantities that provide important information when 
studying confinement and dynamical chiral symmetry breaking in Quantum Chromodynamics 
(QCD). They are necessary inputs for studying bound states and amongst other things, can be 
used to obtain non-perturbative predictions for the running coupling a s . A number of early stud- 
ies [1, 2, 3, 4, 5] found various singularities in the vanishing p 2 limit for the non-perturbative 
propagator dressings. However, more recent results from a variety of methods have found the 
dressings to be finite [6, 7, 8, 9, 10, 11, 12]. 

In a recent study [13], the gluon and ghost propagator dressings, in the absence of quarks, 
were studied. The SDEs were solved self-consistently using a range of approximations for the 
vertices. Emphasis was given to how well these quantities are determined - particularly in the 
physical region around 0. 1 and 1 GeV. It turns out that the vertices are indeed very important if 
numerically precise results are desired. 

Landau gauge is used, since for this kind of problem it is the most appealing theoretically. 
The non-perturbative corrections to the ghost-gluon vertex are simplest in Landau gauge and state- 
ments regarding confinement are typically formulated in Landau gauge also. Furthermore, there 
are several lattice results to which we may compare. 

The dressing functions that we solve for are related directly to the propagators, which for the 
gluon reads, 



where D ilv {p) is the full gluon propagator in Landau gauge and D(p) is the ghost propagator. The 
dressing functions &£ and 9?h contain all of the non-perturbative physics of these two Green's 
functions. The simple transverse dressing of the gluon in Eq. (1.1) is due to the Slavnov-Taylor 
identity for the propagator and is an important feature of the gauge invariance of the theory. 

First we present the solutions of the ghost equation that are obtainable using a fixed gluon 
input. These solutions then provide a natural starting point for investigating the coupled system 
where simultaneous solutions of both propagator dressings are found. During the course of this 
work it became apparent that many simple vertices do not admit self-consistent solutions. We test a 
range of vertices with a range of solutions. A preferred system motivated by theoretical arguments 
is then selected. We conclude by comparing to the extant lattice data. 

2. The Ghost Equation 

Using a fixed gluon input we may solve the Ghost SDE alone and investigate its sensitivity to 
a range of input vertices. This is useful because this type of SDE is very simple to solve, and it 
teaches us what to expect when solving the more complicated, coupled equations self-consistently. 
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In particular, the choice of ghost-gluon vertex applied here does have an effect, and understanding 
this is useful when solving the coupled system of gluons and ghosts. 

In order to study the ghost equation in isolation we are required to provide a gluon input, we 
use the following model [14], 

&e(p 2 ) = tt • (2. i) 

- 2 + ^(l + I^¥Log(^)) E 

In this form the gluon contains a mass term m, which is typically ^(Aqcd) in lattice and other SDE 
studies. The remaining terms arise in order to reproduce the one-loop behaviour. This is similar 
to many contemporary representations of the gluon dressing function, and a closely related form 
has recently been applied to a phenomenological study of hadron physics observables [15, 16]. 
This form reproduces the leading resummed logarithm one finds from perturbation theory at large 
momenta. 
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Figure 1: The ghost-gluon vertex indicating the momentum definition we adopt, the outgoing ghost mo- 
mentum is q and the gluon momentum is k. 

The next step is to specify the ghost-gluon vertex which in principle is all that is required in 
order to solve the ghost equation given a fixed gluon input. The full ghost-gluon vertex has the 
form, 

^ = igf bc {a^p^+Pik^q)^) (2.2) 

where a(k,p,q) and fi(k,p,q) are the non-perturbative dressing functions and the momenta are 
defined in Fig. 1. Initially we make use of the bare vertex in place of the fully dressed one, where 
a = 1 and /3 = 0. This is motivated by Taylor's theorem [17] which constrains the sum of the 
functions to unity at vanishing incoming ghost momentum: p = 0. When solving for the gluon, the 
following extended vertices will be used, 

= igf abc (q^-k^^k^q)) , (2.3) 

if ) = igf abc (qt-p^Jk&R(k,p,q)\ , (2.4) 

where J^ir = 1 for small momenta and J^ir = for large momenta. .Vt\r is a non-perturbative 
normalisation parameter. 

The ghost-gluon vertex in Eq. (2.3) is defined to be transverse in the small-fc limit, and has 
been used to solve these equations previously [8]. In the ghost equation, only the q^ term survives 
and hence this gives identical results to the bare vertex when solving for the ghost alone. This is not 
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the case for , which is a minor modification that has been chosen in order to satisfy the Taylor 
condition. The additional term now vanishes when p — > 0, but contributes in both equations. 

The final requirement is then to specify the renormalisation procedure. In the numerical 
method adopted here, it is most straightforward to use the momentum subtraction scheme, where 
the equations are subtracted from themselves at an arbitrary value of the propagator momentum. 
There are two choices that influence which solutions may be obtained. A subtraction at p 2 = may 
produce a singular ghost dressing function. However a subtraction at some perturbative momentum 
can be applied and in some ways is better since the dressings are already well known in that region. 
The Schwinger-Dyson equations tell us how these evolve from the known results at large momenta 
into the non-perturbative region. The alternative is to make an assumption about the ghost dressing 
at zero momentum. 

The renormalised ghost equation with a cutoff regularisation reads, 

&h-\p 2 )=Z 3 (ii 2 ,A 2 )+n gc (p 2 ,ii 2 ) (2-5) 
= Wh(n 2 ) + U gc (p 2 ^ 2 )+U gc {^ 2 ^ 2 ) (2.6) 

in the second line we eliminate the renormalisation constant Z^jJL 2 ^ A 2 ). The function Yl gc (p 2 ,iJ, 2 ) 
is the ghost self-energy diagram, the details of which are given in [13]. This is the equation that is 
solved numerically. To check the dependence and range of solutions available, we vary the value of 
1 /@h(0). The solutions are given in Fig. 2. We observe two classes of solutions. With Wh(0) > 2 
then the Ghost dressings in the perturbative region are practically indistinguishable. In the small p 2 
region any value in the range oo > £f /j(0) > 2 produces the same perturbative solution. Reducing 
@h(0) there is then some critical value where the solutions in the perturbative region also begin to 
fall. The precise value of h(0) for which this occurs is dependent upon the couplings, the vertex 
and the gluon. All of the solutions are monotonic and decrease with increasing p 2 . By specifying 
the value at p 2 = then all other points are implicitly given by the equation. 
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Figure 2: Examples of ghost solutions on a log-log plot subtracting at zero momentum. The solid curves 
display different ghost dressings @h(p 2 ), while the dashed curve is the gluon dressing ( Sl{jp 2 ). Only the 
specified subnaction value 5f/i(0) is varied between the solutions which can be read from the plot. The 
^h(p 2 = 0) value is fixed and the ghost equation solved with the depicted gluon until the ghost inputs and 
outputs are self-consistent. The units of p 2 are arbitrary since we have not fixed the coupling to the physical 
value, but may be considered to be &(\ GeV 2 ). 
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We may then switch to the other viewpoint of subtracting in the perturbative region and choose 
the typical perturbative condition ^h{pL 2 ) = 1. In doing this, none of the solutions subtracted at 
zero in Fig. 2 are obtained. The ghost dressing rises from unity at p 1 = ju 2 to a finite value at 
zero momentum. It is possible to subtract at zero and fine tune to this value and obtain the same 
solution, although this is a special case. 

Specifying ^/i(/x 2 ) = 1 as one would do perturbatively leads to the interesting consequence 
that only one of these many Ghost dressings is obtained. This is distinct from singular solution and 
importantly no renormalisation group running connects this solution to the singular one. Since the 
expectation is that asymptotically free QCD at large momentum transfers is accurately described 
by the perturbative solution, this is the solution we favour. This effect is depicted in Fig. 3, where 
the solid curve is the unique solution found by imposing ^/j(/i 2 ) = 1 whilst the remaining curves 
are those subtracted at zero from Fig. 2. 

Hence we carry forward this arrangement when solving the gluon equation. It is worth noting 
that this same effect has been tested for gluons that vanish more or less quickly and the same 
qualitative features arise. Similar effects have been observed in other studies [18, 19, 20, 21]. 
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p- 

Figure 3: The dotted curves correspond to the zero momentum subtracted solutions and the colours match 
Fig. 5. The dashed curve is the gluon input, &£(p 2 ), and the solid curve is the physically relevant solution 
of the ghost equation, @h(p 2 ). The units of p 2 are arbitrary. 

Finally we compare the result using which is identical to the bare vertex in the ghost 

(2) 

equation, to the solution found using ; . This comparison is shown in Fig. 4 where a large 
difference is apparent between the two results in the non-perturbative region. 

3. The coupled Gluon and Ghost system 

The coupled system is naturally more complicated, although we already have most of the 
ingredients that are needed. We choose to work using only the one-loop graphs since the two-loop 
graphs are a considerable complication numerically. It is possible that these terms are important 
and in order to obtain numerically precise results it is likely that they are essential. The gluon 
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2.5 r 
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Figure 4: A comparison of two solutions normalised using the perturbative condition. The (red) dotted 
curve is the ghost solution, Gh(p 2 ), with a bare vertex or transverse vertex. The (blue) dashed curve is the 
ghost dressing obtained using The (black) solid curve is the fixed gluon dressing input &l(p 2 ) from 
the model. The units of p 2 are arbitrary internal units. 



equation has the form 

#rV) =z 3 (At 2 ,A 2 )+n 2c (/, A t 2 )+n 2 ,(p 2 , i u 2 ) O-i) 

= ^r\^) + U 2c {p\^)-U 2c {^\^)+U 2g {p\^)-U 2g {^^ 2 ) (3.2) 

where in the second line we have subtracted to remove the renormalisation constant, Yl 2g is the 
gluon loop contribution to the vacuum polarisation, and Yl 2c arises due to the ghost loop. The 
ghost-gluon vertices we have already specified are sufficient for Yl 2c . The final requirement is to 
specify the triple-gluon vertex. To illustrate the uncertainty we select three different vertices and 
check the sensitivity in the solutions. 

The first we apply is a simple form that reproduces the resummed leading logarithmic running 
found at one-loop order in perturbation theory [4], 

r (A) (kpq) _ gft^W) r (0) (k , 3 3 

this is proportional to the bare Lorentz structure F^l p (k,p,q) and p and q are chosen to be the 
internally contracted legs under the loop integration. What one immediately notices about this 
form is that it is not symmetric under interchange of legs. Since this vertex is expected to be 
symmetric to all orders and since the Ward-Slavnov-Taylor Identity (WSTI) is also symmetric 
under the interchange of legs, a symmetric form is desirable. Hence we also apply, 

where the ghost-over-gluon form is suggested by the WSTI, and for simplicity this is chosen to be 
proportional to the bare Lorentz structure. The final dressing of the triple-gluon vertex we apply is 
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a solution of the WSTI using a bare ghost-gluon scattering kernel, 

AC),, 1 ( &h(q 2 ) &h(q 2 ) \ s 

( \ 



1 (&h(k 2 ) <£h(k 2 )\ 



2 \&t(q 2 ) ' 
1 f&h(p 2 ) &h(p 2 ) 



+ ■ 



2 V &l(k 2 ) &£(q 2 



(3.5) 



now each term of the bare structure receives its own non-perturbative dressing and these have a 
symmetric structure. More elaborate forms are also possible and eventually something more will 
be required in order to obtain a triple-gluon vertex that has a consistent ghost-gluon scattering 
kernel with the ghost-gluon vertex as prescribed by the WSTI. 



2.5 r 



Ghost Dressings 




Figure 5: The range of solutions obtained using the possible vertex combinations. The label (i,j) refers to 
the vertices used in obtaining the solutions corresponding to for the ghost-gluon vertex and rj/y p for the 

triple-gluon vertex. The missing curve corresponds to the and T^vp> self-consistent solutions were not 
obtained there. The parameters are not varied between the solutions, only the vertices. The units of p 2 are 
arbitrary since we have not fixed the coupling to a physical value. 

The next step is to implement these vertices numerically in the ghost and gluon SDEs and 
to solve them self-consistently. This is done numerically in the momentum subtraction scheme, 
subtracting both equations at the renormalisation point jU. The equations are solved iteratively 
and are fairly insensitive to the initial function. The model gluon in Eq. (2.1) may be used and 
&h(p 2 ) = 1 is adequate for the ghost. Convergence is achieved rapidly using a Newton-Raphson 
method, although a natural iterative procedure is sufficient. 

In Fig. 5 we see a broad range of gluon dressings given the different vertex inputs. The 
position of the peak of the gluon propagator moves around considerably. This is important because 
it is related to the effective gluon mass and physical quantities depend on this, as has been shown 
in [15, 16]. This plot tells us that just using any old vertex is insufficient and one should appeal to 
all possible constraints when choosing this input. 

These differences in the propagator dressings arise from differences in the vacuum polarisation 
integrals, which come directly from the vertices. It is informative to look inside the contributions 
for the gluon from each loop to disentangle the range of vertex contributions. The simplest triple- 
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Figure 6: Gluon polarization functions, p 2 Yl2 C {p 2 ) (positive values) and p 2 Yl2 g (p 2 ) (negative values). Left: 
IR region, Right: UV region. In the UV the curves are indistinguishable on this scale. In the IR, the 
ghost-loop curves for multiple solutions lie on top of each other. The bare ghost-gluon vertex, not shown 
here, always gives a vanishing IR contribution given an IR vanishing gluon dressing and an IR finite ghost 
dressing. The functions p 2 H2 g {p 2 , fl 2 ) are labelled according to the input vertices with the contributions 
from rL and rl/yp labelled as (i,j) in the plot. The units of p 2 are arbitrary. 



In order to see this we take the subtracted contribution from each loop integration and since they 
typically behave as l/p 2 for small p 2 , we multiply them by p 2 in order to get a better view of the 
differences. This we show in Fig. 6, where the quantities plotted are p 2 (n,(p 2 ,/^ 2 ) — n,(/i 2 ,jU 2 )), 
from Eq. 3.2. In the right panel we see the large p 2 region, which is essentially determined by 
perturbation theory, where all contributions are the same. Differences arise in the left panel of 
Fig. 6 where we zoom in on the small p 2 region. 

The negative curves are the contributions from the triple-gluon vertices. The plotted functions 
p 2 Yl2 g {p 2 , IX 2 ) are labelled according to the input vertices with the contributions from rL and rj/(L 
labelled as (/, j) in the plot. It is evident that T^p leads to a vanishing contribution when plotted 
in this way whilst the sensibly dressed symmetric vertices do not, they are finite and negative. The 
gluon dressing function is not expected to change sign and is positive in the perturbative region, 
so it is necessary that a positive contribution also arises to cancel that from the gluon loop. In the 
one-loop-only system with these triple-gluon vertices, a positive contribution must arise from the 
ghost-loop terms. This is indeed the case since we have self-consistent solutions for the curves that 
are plotted. The combination (1 ,B) is missing since for this set of parameters it did not satisfy the 
condition. It is then not possible to find self-consistent solutions. The same condition is also why 
we need non-trivial ghost-gluon vertices, since the bare vertex gives a vanishing contribution when 
plotted in this way. This leads directly to a change of sign in the gluon dressing and no stable set 
of solutions can be found. 

The fact that a symmetric triple-gluon dressing results in a negative contribution is very im- 
portant, since it leads to the condition that the ghost loop, and hence the ghost-gluon vertex, need 
not be transverse alone and can work together with the other loops, as they do in the perturbative 
region, in order to produce a transverse gluon dressing. 
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4. Comparison with Lattice QCD 

Primarily we have been motivated by theoretical issues encountered in solving the Schwinger- 
Dyson equations for the gluon and ghost propagators. However lattice studies exist for these quan- 
tities and the results are complementary. The possible issues that result from a applying a finite 
grid to Euclidean spacetime are quite different to those that we may induce here by truncation, so 
a comparison is a useful independent cross-check. Importantly there are lattice computations that 
are precise and in the pure gauge sector. 

An early lattice result favoured a finite, massive solution of the gluon [22]. Following this 
many recent lattice studies now exist using very large lattices [6, 11, 12, 23]. We compare our 
calculations to [9], which provides results in Landau gauge for both dressing functions. The quali- 
tative behaviour is the same, with a finite ghost dressing function and a vanishing gluon propagator 
dressing function. 



12r 




0.001 0.1 10 1000 h)- 5 0.001 0.1 10 1000 

p p 



Figure 7: Self-consistent solutions tuned to lattice solutions, showing the dressing functions. The solid 
curves depict a smooth fit to the lattice data. The heavier region is where the functions are represented 
by lattice data and the feint region represents the natural extrapolation. The broken curves are the tuned 
solutions. Left: The blue curve corresponds to the gluon dressing function &£(p 2 ), this vanishes as p 2 —> 0, 
the red curve corresponds to the ghost dressing function &h(p 2 ). Right: The blue curve corresponds to 
&e(p 2 ) jp 2 , this is ~ 10 as p 2 0, the red curve corresponds to the ghost dressing function again. 

In Fig. 7 we show the solutions of the equations obtained self-consistently with three different 

sets of parameters that reproduce different features of the lattice data. In obtaining these solu- 

(c) 

tions we use the preferred set of vertices, r^ vp from Eq. (3.5) and the ghost-gluon vertex given in 
Eq. (2.4). We then tune the parameters to obtain a reasonable representation of the lattice data. We 
are not able to find a close fit over the whole momentum region. All three sets of parameters fail 
to reproduce the magnitude of the peak of the gluon propagator dressing function shown in the left 
panel of Fig. 7. It is possible that this is due to the effects of the two-loop graphs. However to be 
sure one would need to perform the full dressed two-loop integrals. 

Another difference is the value of the coupling in the perturbative region. In all three sets of 
solutions, the lattice appears to have a larger value of the coupling. This leads to a slower increase 
in the SDE solutions between the peak of the gluon propagator and the subtraction point. A similar 
effect is seen in the ghost dressings in this region. In the non-perturbative region it is possible to 
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find similar values of the dressings simultaneously: The dashed solution is the closest we show and 
the effect is most visible in the right panel of Fig. 7. 

5. Conclusion 

We solved the Schwinger-Dyson Equations of QCD in the Landau gauge in the absence of 
quarks. In solving the ghost equation we found the singular ghost is excluded when a subtraction 
point in the perturbative region is selected. This was true for all vertices that we tested and for a 
range of gluon inputs although only one was shown. This agrees with earlier findings [18]. 
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Figure 8: The running coupling aj(p 2 ) is shown as the black solid curve. The running coupling (Xec{p 2 ) 
is shown as the green solid curve. The gluon dressing %?l(p 2 ) is the dashed blue curve and the dotted red 
curve is the ghost dressing function @h(p 2 ). The units of p 2 are now very close to GeV 2 . 

When solving the gluon equation we found that non-trivial vertex dressings were essential in 
that simple vertices did not lead to self-consistent solutions. It was found that symmetric triple- 
gluon vertex dressings lead generically to a non-zero contribution of the gluon-loop diagram in the 
non-perturbative region. The main implication of this is that the ghost-loop term is not required to 
be transverse alone and that both terms can contribute to the transversality of the gluon propagator, 
as they do in the perturbative region. This combined with Taylor's theorem led us to test a simple 
modification of a vertex already known to give self-consistent solutions with in Eq. (2.4). 
This modification protects Taylor's theorem whilst sacrificing the transversality of that individual 
loop. It is interesting to note that another ghost-gluon vertex dressing presented at this workshop 
leads to similar effects in the ghost equation [24]. 

A range of solutions were found by selecting a conservative set of parameters and testing the 
effects of using the different vertices. A wide variation was found in the physical region leading us 
to the suggestion that more advanced vertices are required in order to obtain numerically precise 
results. Proposals for new vertices in a related method specifically investigating the mass generation 
mechanism have been suggested [25, 26]. 

We may use the method we present here to produce a running coupling by considering the 
ghost-gluon vertex renormalisation. There are two definitions we apply, first using the standard 
Taylor definition which we denote aj and a second form, CCec related to the Taylor form, but 
removes the mass term from the gluon dressing [27]. These curves and a set of solutions fixed 
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using physically meaningful parameters are given in Fig. 8. The effective gluon mass we found 
from this gluon propagator was m g = 400 MeV. 

These equations produced self-consistent solutions of the ghost and gluon propagator dress- 
ings that were in qualitative agreement with lattice QCD results. Quantitative differences from this 
method could arise in either the vertex dressings or two-loop graphs. Regardless of precision, qual- 
itatively all results are the same and they point to a finite gluon propagator in the non-perturbative 
region, corresponding to a dynamically generated gluon mass as proposed by Cornwall [28]. 
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